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Abstract 


In natural SUSY models higgsinos are always light because jj? cannot be much larger 
than M^, while squarks and gluinos may be very heavy. Unless gluinos are discovered 
at LHC13, the commonly assumed unihcation of gaugino mass parameters will imply 
correspondingly heavy winos and binos, resulting in a higgsino-like LSP and small inter- 
higgsino mass splittings. The small visible energy release in higgsino decays makes their 
pair production difficult to detect at the LHC. Relaxing gaugino mass universality allows 
for relatively light winos and binos without violating LHC gluino mass bounds and without 
affecting naturalness. In the case where the bino mass Mi ~ g, then one obtains a 
mixed bino-higgsino LSP with instead sizable Wi — Zi and Z 2 — Z\ mass gaps. The 
thermal neutralino abundance can match the measured dark matter density in contrast 
to models with a higgsino-like LSP where WIMPs (weakly interacting massive particles) 
are underproduced by factors of 10-15. If instead M 2 ~ g, then one obtains a mixed 
wino-higgsino LSP with large Z 2 — Z\ but small W\ — Z\ mass gaps with still an under¬ 
abundance of thermally-produced WIMPs. Portions of the light wino parameter space 
may already be excluded by the IceCube upper bound on the spin-dependent neutralino- 
nucleon cross section. We discuss dark matter detection in other direct and indirect 
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detection experiments and caution that the bounds from these must be interpreted with 
care. Finally, we show that LHC 13 experiments should be able to probe these non- 
universal mass scenarios via a variety of channels including multi-lepton -|- events, 

WZ ^ events, Wh + Fi™®® events and -|- Fi™®® events from electroweak 

chargino and neutralino production. 
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1 Introduction 


Results from the first extended runs of LHC at y/s = 7 and 8 TeV have led some authors to 
imply that there is a crisis in supersymmetry (SUSY) phenomenology [T|: how can it be that 
the Higgs and vector boson masses - whose values are related to weak scale soft SUSY breaking 
(SSB) parameters and to the superpotential parameter /i - are clustered near 100 GeV while 
superpartner masses, whose values are also determined by soft SUSY breaking terms, are so 
heavy that they are beyond the reach of LHC? The superpotential higgsino mass parameter /i 
and the SSB Higgs mass parameters enter via the tree-level Higgs potential, whereas other SSB 
parameters - specihcally, those that affect sparticles with the largest couplings to the Higgs 
sector - only enter at higher order. This is clearly evident, for example, in the well-known 
expression, 

Mz _ + ^d- tan^ /? 2 ^ 

2 tan2/?-l ^ 

for the Z mass, where and denote the 1-loop corrections explicitly given in the Appendix 
of Ref. [2j. SUSY models requiring large cancellations between the various terms on the right- 
hand-side of ([^ to reproduce the measured value of M| are regarded as unnatural, or hne- 
tunedfl 

Several measures have been proposed [31 n El El E] to quantify the degree of hne-tuning. A 
common feature of these is that they all regard the model to be hne-tuned if /i^ > M|. This is 
because in most models /i directly enters Eq. ([^ as an independent parameter, and unexplained 
cancellations have then to be invoked to obtain the observed value of Mz- In contrast, in most 
models the SSB masses are obtained in terms of a one (or more) model-parameters and so 
are not independent, allowing for the possibility of large cancellations that is ignored by the 
commonly used large log measure. It is a neglect of these parameter-correlations that has led 
some authors to conclude that light top-squarks are a necessary feature of natural SUSY. In fact, 
as we have just argued, it is |/i| ~ Mz and concomitantly the existence of light higgsino^ (and 
not light stops) that is the robust conclusion of naturalness considerations. The importance of 
low pi for electroweak naturalness was recognized by Chan, Chattopadhyay and Nath [S] over 
hfteen years ago and has recently been emphasized in Refs. and by Martin [10] . 

In earlier papers we have developed the radiatively-driven natural SUSY (RNS) framework 
characterized by values of the parameter Aew = 10 — 30 range corresponding to 3-10% elec¬ 
troweak hne-tuning ilElEI. Within this framework, 

• the superpotential pi term has magnitude |/i| ~ 100 — 300 GeV (the closer to Mz the 
better); 

• the up-Higgs soft term is driven radiatively to small negative values m'jj {weak) ~ 

^We emphasize that for superpartners up to a few TeV range, the degree of fine-tuning that we are talking 
about is many orders of magnitude smaller than in the Standard Model because scalar masses do not exhibit 
quadratic sensitivity to physics at the ultra-high scale if SUSY is softly broken. 

^We assume here that there is no SUSY-breaking higgsino mass term (such a term would lead to hard SUSY 
breaking - and so would be automatically forbidden ~ if higgsinos had superpotential Yukawa couplings to any 
Standard Model singlet 0) so that the higgsino mass comes only from the superpotential parameter p,- This is 
the case in all models that we know of. 
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• the magnitude of radiative corrections contained in should be smaller than or com¬ 
parable to M|. This latter condition occurs for TeV-scale highly mixed top squarks- a 
situation which also lifts into the 125 GeV regime [ 6 ]. In contrast, the terms and 

can occur at the multi-TeV level since they are suppressed by tan^/3 where tan/? is 
required to be in the 3-50 range. 

• Since the gluino mass feeds into the stop masses via RG evolution- and thus into 

then low Aew also requires an upper bound on nig A 4 — 5 TeV |2]. Of course, M 3 is also 
bounded from below by the experimental bound of nig ~ 1.3 TeV based on LHG 8 searches 
within the context of SUSY models like mSUGRA/GMSSM [11] or within simplihed 
models [T2|. 

• First and second generation sfermions can be allowed anywhere in the ~ 5 —20 TeV range 
without jeapordizing naturalness [13]. The higher range of values ameliorates the SUSY 
flavour, CP, gravitino and proton-decay problems due to decoupling. 

Inspired by gauge coupling unihcation, in these previous studies we had assumed gaugino mass 
unihcation as well as naturalness. From gaugino mass unihcation one expects at the weak scale 
that Ml ~ M 3/7 and M 2 ~ 2 M 3/7 so that the LHG 8 lower bound on M 3 also provides a lower 
bound on Mi and M 2 . In this case, for natural SUSY which respects LHG 8 bounds, we expect 
the mass hierarchy |/r| < Mi < M 2 < M 3 to occur. Thus, in the RNS model which we take as 
the paradigm case for the study of natural SUSY, one expects four light higgsino states with 
mass U 7 .^±, where the lightest higgsino Zi acts as the lightest-SUSY-particle or 

LSP. In particular, mixed higgsino-bino or higgino-wino LSPs are not allowed if the gluino is 
heavy. 

Gollider signals as well as cosmology depend sensitively on the nature of the LSP. For 
instance, in the RNS framework with gaugino masses near the TeV range, we expect the light 
electroweak -inos and Zi ^2 to be dominantly higgsino-like with typically small 
and mg^ — mg^ mass splittings of order 10-20 GeV [2]. Such a small mass splitting results in 
only soft visible energy release from the heavier higgsino three-body decays to the Zi. This 
situation makes pair production of higgsinos very difficult to detect at LHG [TH UHl UHl 113 [TH] 
in spite of their relatively small masses and correspondingly large production cross sections; 
other superpartners may be very heavy, and possibly beyond the reach of the LHG. In contrast, 
in models with light gauginos and heavy higgsinos, the mass gap between the bino and wino- 
like states tends to be large (if gaugino mass unihcation is assumed), and signals from wino 
pair production followed by their decays to bino-like LSPs should be readily detectable. The 
celebrated clean trilepton signature arising from W 1 Z 2 production is perhaps the best-known 
example. 

The phenomenology of dark matter is even much more sensitive to the content of the LSP. 
Higgsino and wino-like LSPs lead to an under-abundance of thermally-produced LSPs whereas 
a bino-like LSP leads to overproduction of WIMPs (weakly interacting massive particles) unless 
the neutralino annihilation rate is dynamically enhanced, e.g. via an s-channel resonance or 
via co-annihilation, or their density is diluted by entropy production late in the history of the 
Universe. In the wino- or higgsino-LSP cases, if one solves the strong CP problem via a quasi- 
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visible axion [T^, then the dark matter is expected to occur as an axion-neutralino admixture, 
i.e. two dark matter particles 1201 . 

Gaugino mass unihcation - well-motivated as it may be - is by no means sacrosanct. Phe¬ 
nomenologically, while the high scale value of M 3 is required to be large by LHC 8 constraints 
on rrtg, Mi and/or M2 may well have much smaller magnitudes without impacting naturalness. 
These considerations motivated us to examine how the phenomenology of natural SUSY models 
with Ip I ~ 100 — 300 GeV may be altered if we give up the gaugino mass unihcation assump¬ 
tion and allow for the possibility that the bino or/and wino also happens to be light. The LSP 
(and possibly also other electroweak-inos) would then be mixtures of higgsinos and electroweak 
gauginos, or may even be very nearly bino- or wino-like, resulting in very different mass and 
mixing patterns from expectations within the RNS framework. A mixed bino-higgsino LSP 
could well lead to the observed relic-density for thermally produced neutralinos. We acknowl¬ 
edge that small values of gaugino mass parameters would have to be regarded as fortuitous 
from the perspective of naturalness. Nevertheless since light winos/binos do not jeapordize 
naturalness, in the absence of any compelling theory of the origin of SSB parameters, we felt a 
phenomenological study of this situation is justihed by our philosophy that it is best to “leave 
no stone unturned” in the search for natural SUSY at the LHG. 

Non-universal gaugino masses (NUGM) can occur in GUT models wherein the gauge kinetic 
function transforms non-trivially as the direct product of two adjoints [211 [ 22 ]. Or, it may be 
that GUTs play no role, and that unihcation occurs within the string-model context. Models 
with mixed anomaly- and gravity-mediation contributions to gauginos masses also lead to non- 
universal gaugino mass parameters [23] • Investigation of how the phenomenology of natural 
SUSY models is modihed from RNS expectations forms the subject of this paper. Naturalness 
in the context of non-universal gaugino masses has also been considered in Refs. [23] and [10] . 

1.1 Natural SUSY benchmark scenarios 

We begin by exhibiting a sample benchmark point within the framework of the canonical 2- 
extra-parameter non-universal Higgs model (NUHM2) with unified gaugino mass parameters 
and a higgsino-like LSP under the column RNSh in Table [^ This point has parameters mo = 
5000 GeV, mi /2 = 700 GeV, Aq = —8000 GeV and tan/3 = 10 with (/i, tua) = (200,1000) 
GeV. The RNSh point has Aew = 9-6 corresponding to about 10% electroweak hne-tuning, 
and mh = 124.3 GeV while rrig ~ 1.8 TeV with m^ = 5.2 TeV. It is safely beyond LHG 8 reach. 
The lightest neutralino is dominantly higgsino-like (higgsino-wino-bino composition is listed as 

and dehned similarly to Ref. [23]) and has mass m^^ = 188 GeV 
and thermally-produced neutralino relic density [23] = 0.013. SUSY contributions to the 

branching fraction for 6 —)■ 57 are negligible so that this is close to its SM value [2Z] and in accord 
with experiment [2H]- The spin-independent neutralino-proton scattering cross section shown in 
the third-last row of the table naively violates the bound a^^(Zip) ~ (2 — 3) x 10“® pb from the 
LUX experiment [22], but we note that this bound is obtained assuming that the neutralino 
comprises all of the cold dark matter. In our case, the thermal neutralino contribution is 
just about 10% of the total DM contribution, and this point is in accord with the constraint 
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upon scaling the expected event rate by ^ /Q.12. Also shown is the spin-dependent 

neutralino-nucleon scattering cross-section which is on the edge of observability of the IceCube 
experiment [3l] which is searching for high energy neutrinos from neutralinos captured by the 
sun and annihilated in the solar core. Since the neutralino density in the sun is obtained 
assuming equilibration between the capture rate and the annihilation rate of neutralinos, no 
scaling is needed in this case. The other columns display natural SUSY benchmark points 
where the bino or the wino mass parameters are dialed to relatively low values resulting in 
natural SUSY models with either a bino-like (RNSb) or wino-like (RNSw) LSP. These cases 
will be discussed in detail in the following sections. 

1.2 Remainder of paper 

The remainder of this paper is organized as follows. In Sec. we hrst investigate the case 
of \Mi\ ~ |/i| M 2,3 where we treat Mi as an additional phenomenological parameterj^ In 

this case, the LSP can become mixed bino-higgsino or even mainly bino-like. Note also that 
while we can always choose one of the gaugino mass parameters to be positive, the signs of 
the remaining ones are physical. In our study, we will examine both signs of the gaugino mass 
parameters that are assumed to depart from universality. In Sec. we investigate the case 
with IM 2 I ~ |/u| Ml ,3 which can generate a wino-like LSP. In Sec. we examine the more 
general case where both |Mi| and IM 2 I are simultaneously comparable to |/i|. Our conclusions 
are presented in Sec. 


2 Natural SUSY with a bino-like LSP 

In this section, we examine how the phenomenology of natural SUSY models is altered if 
we allow for non-universal gaugino mass parameters, and let the GUT scale bino mass vary 
independently. To this end, we adopt the RNSh benchmark point from Table [TJ but now allow 
Ml to be a free parameter, positive or negative. To generate spectra and the value of Aew, 
we adopt the Isajet 7.84 spectrum generator |32]. In Fig. we show by red circles the value 
of Aew versus the GUT scale value of Mi. We see that - aside from numerical instabilities 
arising from our iterative solution to the SUSY RGBs - the value of Aew stays nearly constant 
so that, as anticipated, varying Mi hardly affects the degree of electro-weak hne-tuning. 

In Fig. 1^ we show the mass values of the charginos and neutralinos as Mi is varied between 
-700 GeV to 700 GeV. For Mi = 700 GeV, the gaugino mass unihcation point, we hnd that 
Wi and Zi ,2 are all higgsino-like with mass values clustered around p = 200 GeV while the 
bino-like Z 3 lies near 300 GeV and the wino-like Z 4 and W 2 lie at ~ 600 GeV. As Mi is lowered, 
then the bino component of Zi increases while the bino-component of Z 3 decreases. The mass 
eigenvalues track the gaugino/higgsino content, and as we pass through Mi = 300 GeV, the Zx 

^We remark that other processes may further alter the neutralino relic density from its thermal value, 
increasing it if there are late decays of heavy particles to neutralinos, or diluting it if these decay into SM 
particles. For more detailed discussion of non-thermally-produced dark matter, see the recent review |30j . 

^We frequently denote both the GUT and weak scale values of the gaugino mass parameters by Mj. We 
assume that it will be clear from the context which case is being used so that this abuse of notation will not 
cause any confusion. 
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parameter 

RNSh 

RNS5 

RNSw 

Ml (GUT) 

700 

380 

700 

M 2 (GUT) 

700 

700 

175 

Ms (GUT) 

700 

700 

700 

mg 

1795.8 

1796.2 

1809.8 

mu^ 

5116.2 

5116.2 

5100.7 

rriu^ 

5273.3 

5271.3 

5277.4 

me« 

4809.0 

4804.4 

4806.7 


1435.1 

1438.1 

1478.3 


3601.2 

3603.3 

3584.9 


3629.4 

3631.5 

3611.6 


5003.9 

5003.6 

5007.4 

rrif^ 

4735.6 

4731.1 

4733.9 

mf2 

5071.9 

5070.8 

5053.9 


5079.2 

5078.1 

5060.8 

^W2 

610.9 

611.0 

248.4 


205.3 

205.3 

121.5 


621.4 

621.5 

322.1 


322.0 

217.9 

237.8 

-zY 

209.3 

209.8 

211.8 


187.8 

149.5 

114.2 

rrih 

124.3 

124.2 

124.3 

VTU 

0.96 

0.57 

0.60 


-0.14 

0.07 

-0.80 


0.24 

-0.82 

0.08 

Aew 

9.6 

9.6 

10.8 

Zi 

0.013 

0.11 

0.0015 

BF{b —)■ s'j) 

3.3 X 10-^ 

3.3 X 10-^ 

3.3 X 10-^ 

a^\Zip) (pb) 

1.6 X 10-® 

1.7 X 10“® 

4.3 X 10-® 

a^^{Zip) (pb) 

1.7 X 10-^ 

2.8 X 10“^ 

8.9 X 10-^ 

(an)|„^o (cmVsec) 

2.0 X 10-25 

1.8 X 10-26 

1.7 X 10-24 


Table 1; Input parameters and masses in GeV units for three Natural SUSY benchmark points 
with /i = 200 GeV and = 1000 GeV. We also take mo = 5000 GeV, Aq = —8000 GeV and 
tan (5 = 10. Also shown are the values of several non-accelerator observables. 
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Figure 1: Variation in fine-tuning measure Aew vs. Mi (red circles) or M 2 (blue pluses), with 
all other parameters fixed at their values for the RNS SUSY benchmark model point in Table [T} 
Here, and in subsequent figures the Mj on the horizontal axis is the value of the corresponding 
gaugino mass parameter renormalized at the GUT scale. We cut the graphs off if the lighter 
chargino mass falls below 100 GeV. 


and Z 3 exchange identities and interchange from being bino-like to higgsino-like. A similar level 
crossing is seen on the negative Mi side of the hgure. Since there is no charged bino, the values 
of ^ remain constant (at /i and M 2 (weak)) with variation of Mi. Since the value of is 
decreasing as Mi decreases, then the mass gaps ^Z 2 ~ ^Zi increase. The 

mass gaps reach values of ~ 150 GeV for Mi as small as 50 GeV. This should render signals 
from W 1 Z 2 and WiWi production much easier to detect at the LHG as compared to the RNSh 
case. 

In Fig. we show the thermally-produced neutralino relic density as calculated using the 
IsaReD program [26]. The value of begins at ~ 0.01 for |Mi| = 700 GeV which is 

typical for a higgsino-like LSP of mass 200 GeV. As |Mi| decreases, then the bino content of 
Zi becomes larger - reducing the annihilation cross section - so that the thermal relic density 
correspondingly increases. For |Mi| ~ 380 GeV, the value of reaches 0.12, i.e. it saturates 

the measured DM abundance, and we have the so-called well-tempered neutralino. For even 
lower values of |Mi|, then neutralinos are unable to annihilate efficiently and exceeds 1 

except for special values where the neutralino annihilation cross-section is resonance-enhanced. 
For |Mi| ~ 150 GeV, then the bino-like neutralino has mass ~ mh/‘^ so that neutralinos 
can efficiently annihilate through the light Higgs resonance. The annihilation rate at resonance 
is not quite symmetric for the two signs of Mi. For even lower values of |Mi|, then ~ M^/2 
so that neutralinos efficiently annihilate through the Z boson pole. At values of |Mi| < 100 
GeV we move below the Z-resonance and due to the increasing bino content of Zi, the LSP 
annihilation cross section becomes even smaller, leading to an even larger thermal relic density]^ 

®We remind the reader that these parameter regions with seemingly too large a thermal neutralino relic 
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Figure 2: Variation of electroweak-ino masses vs. Mi for a general RNS SUSY benchmark 
model with variable Mi and M 2 = M 3 


We display the SUSY spectrum for Mi (GUT) = 380 GeV, the value for which the thermal 
neutralino relic density essentially saturates the measured abundance so that = 

0.12, in Table as RNSb. In this case, the Zi is a bino-higgsino admixture albeit already it is 
dominantly bino-like. The mass gap is ~ 56 GeV while the mass gap is 

~ 60 GeV. 

2.1 Implications for LHC13 

The possibility of non-universal gaugino mass parameters has important implications for dis¬ 
covery of natural SUSY at LHG13. 

2.1.1 Gluino pair production: multi-jet -|- U™®® events 

Since squarks are very heavy, the multijet -|- U™'®® signal mainly arises from pp —)■ ggX followed 
by gluino cascade decays mainly via g —)• tbWj and tiZi. For a hxed rrig, but varying Mi, one 
still expects multi-lepton plus multi-jet-|-U™®® events at a rate which mainly depends on the 
value of rrig. For discovery via gluino pair production, the LHG13 reach - which extends to 
about rrig ~ 1.7 TeV (for nig -C rriq) for 100 fb“^ of integrated luminosity |T1] - tends to be 
dominated by multi-jet-t-U™^®® channel and so changes little compared to the case of universal 
gaugino masses. For the RNS point in question, the gluino dominantly decays via g —)■ Uf, and 
the ti subsequently decays via ti —)■ bWi,tZi^ 2 , 3 - Within the gluino pair cascade decay events, 
the isolated multi-lepton content should increase with decreasing Mi due to the increased 

density should not summarily be excluded because the neutralino relic density can be diluted if, for instance, 
there are heavy particles with late decays into Standard Model particles in the early universe. 
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Figure 3: Variation of vs. Mi for a general RNS SUSY benchmark model with variable 

Ml and M 2 = M 3 . The dashed line shows the measured value of the cold dark matter relic 
density. 

mass gap between Wi — Zi and ^ 2,3 — Zi since one may also obtain energetic leptons from 
Wi —>■ ^l'(,Zl and Z 2 —)■ Zl^'^^~ three body decays in addition to those from top or Z 3 decays. 
If Ml is sufficiently small, then the two-body decays Wi ZiW and Z 2 —)■ ZiZ, Zih open 
up. The latter two decays, if open, tend to occur at comparable rates in natural SUSY with 
a bino-like LSP since the lighter -inos tend to be a gaugino-higgsino admixture. The isolated 
opposite-sign/same flavor (OS/SF) dileptons present in cascade decay events will have mass 
edges located at for three-body decays, or else real Z —)■ or h —)■ 66 pairs will 

appear in the case of two-body decays of Z 2 and Z^\ 

2.1.2 Electroweak -ino pair production 

For electroweak-ino pair production, allowing non-universality in the gaugino sector changes the 
situation quite dramatically. In the case of RNS with gaugino mass unihcation, the higgsino pair 
production reactions pp —)■ IVi^Wf and W 1 Z 12 are largely invisible due to the small mass gaps 
la. It may, however, be possible to detect higgsino pair production making use of initial state 
QCD radiation and specially designed analyses if the higgsino mass is below ~ 170 — 200 GeV, 
depending on the integrated luminosity |T 3 IT 8 |^ 

The wino pair production process pp —)■ W 2 Z 1 X can lead to a characteristic same-sign 
diboson signature | 33 ] arising from —)■ ZiW^ and Z 4 —)■ decays, where the higgsinos 

decay to only soft visible energy and are largely invisible. 

In contrast, as Mi diminishes, then the growing lUi — Zi and Z 2 — Zi mass gaps give rise 
increasingly to visible decay products and a richer set of electroweak -ino signals. In Fig. 
we show the NLO cross sections obtained using Prospino [Mj for various electroweak-ino pair 
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production reactions versus variable MiiGUT) for the RNS benchmark casej^ As Mi falls to 
lower values, the chargino pair rates remain constant since fi and M 2 do not change. The W 1 W 2 
cross section in the topmost frame is small because squarks are very heavy, and the ZW 1 W 2 
coupling is dynamically suppressed. Although the Wi —)■ ffZi decay products become more 
energetic with reducing |Mi|, the chargino pair signals are typically challenging to extract from 
large SM backgrounds such as W~^W~ production. 

For WiZi ^2 production, the cross sections can be large but the decays give only soft visible 
energy for Mi ~ 700 GeV. But as Mi is lowered, the cross section for W 1 Z 2 remains large but 
the mass gaps increase. Ultimately, the clean trilepton signature should become visible against 
SM backgrounds |3Sl EH]- Also, the reaction pp —)■ W 1 Z 3 has an increasing cross section as 
Ml decreases and should give rise to i + Z events: trileptons where one pair reconstructs a 
real Z 1371, as is the case for the RNSb benchmark point: see also Ref’s. [3H1 EH]. Ultimately, 
the Z 3 —)■ Zih mode also opens up, reducing the trilepton signal but potentially offering an 
opportunity for a search via the Wh channel jlO]. 

In models with heavy squarks, higgsino pair production reactions make the main contribu¬ 
tion to neutralino pair production processes. In many models, |p| is large, making neutralino 
pair production difficult to see at hadron colliders. Natural SUSY models with non-universal 
gaugino masses are an exception as can be seen from the bottom frame of Fig. where we 
show cross-sections for various neutralino pair production processes. The bino-higgsino level 
crossing that we mentioned earlier is also evident: for large Mi the Zi and Z 2 are higgsino-like 
states and Z 1 Z 2 production (solid squares) dominates, whereas for small Mi then Z 2 and Z 3 are 
higgsino-like and Z 2 Z 3 production (left-pointing triangles) is dominant even though the Z 1 Z 2 
and Z 1 Z 3 reactions are kinematically favoured. Also Z 1 Z 2 and Z 2 Z 3 production can lead to 
dilepton and four-lepton hnal states which may be visible, and to ZZ, Zh and hh -|- hnal 
states if \Mi\ is sufficiently small. 

2.2 Implications for ILC physics 

The prospects for SUSY discovery and precision measurements in the RNS model have been 
examined for an International Linear e’''e“ Collider (ILC) with y/s ~ 250 — 1000 GeV in Ref. 
|41j . Such a machine is a higgsino factory in addition to a Higgs factory and even with small (10 
GeV) inter-higgsino mass gaps, SUSY signals should stand out above SM backgrounds. The 
clean environment, together with the availability of polarized electron beams, also allows for 
precision measurements that point to the higgsino origin of these events. The main reactions 
of import are —)■ and Z 1 Z 2 production. 

In the case where Mi is low enough so that one obtains a bino-like LSP, the second higgsino 
state Z 3 also becomes accessible, and reactions involving Z^ provide even richer prospects for 
SUSY discovery. Various SUSY pair production cross sections are shown in Fig. [^versus vari¬ 
able Ml and for y/s = 500 GeV. The electron and positron beams are taken to be unpolarized 
in this hgure. Once again the level crossings between bino and higgsino-like states are evident. 
For the case of unihed gaugino masses with Mi = 700 GeV, then indeed only WiWi and Z 1 Z 2 
are available. However, as \Mi \ is lowered, then ct(IUiIUi) remains constant although the decay 

®Since, as we saw in the previous figures the mixing patterns are roughly symmetric about Mi = 0, and 
because it is relatively time-consuming to run Prospino, we show results only for positive values of Mi. 
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Figure 4: Electroweak -ino pair production cross sections versus Mi for the RNS SUSY bench¬ 
mark model with variable Mi but with M 2 = M 3 
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Figure 5: Chargino and neutralino production cross sections at a linear e’''e“ collider with 
^/s = 500 GeV with unpolarized beams for the RNS SUSY benchmark model with variable Mi 
but with M 2 = M 3 

products of Wi become more energetic once the LSP becomes bino-like and lighter than the 
higgsino. The dijet mass spectrum from Wi —)■ Ziqq' decay allow for precision extraction of 
and and also extraction of the weak scale SUSY parameters /i and also Mi, if the 
bino mass is small enough ga sal E]. 

Turning to neutralino production, we see that higgsino pair production - Z 1 Z 2 production 
if |Mi| is large, and ^ 2-^3 production for small values of |Mi| - dominates the neutralino cross 
section just as in the LHC case. Notice that for 0 < Mi < 300 GeV, .^ 1.^3 production also 
occurs at an observable rate, falling with reducing Mi because of the increasing bino content 
of ZiQ We have checked that the strong dip in a{ZiZ^) around Mi ~ 500 GeV is due to 
an accidental cancellation in the ZZiZ^ coupling]^ -^ 2^3 and Z 1 Z 3 production should lead to 
interesting event topologies, including Z + and h + events where the missing mass 
does not reconstruct to M^, depending on the decay modes of the neutralinos. On the negative 
Ml side, the ZiZ^ cross section is small, except beyond the level crossing at Mi ~ —600 GeV. 

Before closing, we note that these neutralino and chargino cross sections are also sensitive 
to beam polarization. This can serve to extract the gaugino/higgsino content of the charginos 
and neutralinos that are being produced. 

^This is somewhat different from the behaviour in Fig. where we see, for example, that a{ZiZ 2 ) increases 
with reducing Mi. We attribute this to the reduction in mass of the .^ 1.^2 system and the concomitant increase 
of the parton densities at the LHC. 

®The alert reader may wonder why there is no similar dip in Fig. We remark that the code used to make 
Fig. [5] uses tree-level masses and mixings among charginos and neutralinos, whereas Fig. includes effects of 
radiative corrections to the spectrum. These corrections, of course, shift the location as well as depth of the 
dip. We have checked that the coupling is indeed suppressed even with radiative corrections, butj^fwre is no big 
dip, at least within resolution of the scan. Since it has no implications physics-wise because < 7 ( 21 ^ 3 ) in Fig. Id 
is already very small for Mi ^ 500 GeV, we have not attempted to refine this figure. 
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2.3 Implications for dark matter seaches 

In the RNS model with unihed gaugino masses and a higgsino-like LSP, the relic density of 
thermally produced neutralinos is much smaller than the observed density of cold dark mat¬ 
ter. This allows for a contribution from axions [19] that must be present if nature adopts the 
Peccei-Quinn solution to the strong CP problem. In the case of DFSZ axions [H], one also 
gains a solution to the SUSY /i problem and can allow for a natural value of /i ~ 100 — 200 
GeV via radiative PQ breaking Ha. In such models, the DM tends to be axion-dominated 
[16] with a local abundance of neutralino WIMPs reduced by factors of 10-15 from usual ex¬ 
pectations. The reduced local abundance makes direct detection more difficult since detection 
rates depend linearly on the local neutralino abundance. Indirect detection rates from WIMP 
halo annihilations depend on the square of the local abundance so are even more suppressed in 
models where the WIMPs only make up a fraction of the dark matter 071. 

For the more general model where |Mi| may be lower than expected from gaugino mass 
unihcation, the thermally-produced neutralino abundance is increased, and consequently one 
expects a greater fraction of neutralino dark matter compared to axions, assuming there are no 
other processes that affect the neutralino relic density. The increased local neutralino abundance 
leads to more favorable prospects for WIMP direct and indirect detection. 

The spin-independent (SI) WIMP-proton scattering cross section from IsaReS jUj is shown 
in Fig. The curve with red dots shows the case of variable Mi. As Mi decreases from large, 
positive values, then the LSP becomes more of a bino-higgsino admixture. Since the SI cross 
sections proceeds mainly through light Higgs h exchange, and the Higgs-neutralino coupling 
is proportional to a product of gaugino times higgsino components 1251. then the SI direct 
detection cross section increases by up to a factor of ~ 2 for lowered Mi. As Mi is lowered 
even further, then the LSP becomes more purely bino-like, and the SI direct detection cross 
section drops sharply. The sharp dip at Mi ~ —110 GeV is due to the reduction of the hZiZi 
coupling, and also the cancellation between the neutralino scattering through the exchange of 
the light GP-even Higgs and that through the exchange of the heavy GP- even Higgs, denoted 
as the blind spot in dark matter direct detection [25] H9] [50] . The kink at Mi ~ —600 GeV 
occurs due to a change in the composition of the LSP: we see from Fig. |^that the levels are 
getting very close, and the -inos may be switching composition. 

The reader may be concerned that the cross-section in Fig. [^seemingly violated the upper 
limits from LUX (Ref. [20]) of ~ (1 — 2) x 10“® pb for neutralinos in the mass range 20-200 GeV. 
As mentioned previously, we should remember that these limits assume that the LSP saturates 
the observed density of cold dark matter, which is certainly not the case for a higgsino-like LSP 
(large |Mi| values in the hgure). Scaling by the expected fraction of the thermal relic density 
makes the large |Mi| region safe, though on the edge of observability of the LUX experiment, 
if thermal production is assumed to be the complete story of the neutralino relic density. For 
smaller values of |Mi|, where it may also appear that the direct detection bound is violated, this 
clearly is not the case. We should, however, keep in mind that for these ranges of Mi, the direct 
detection rate from which the bound in Ref. [29] is inferred cannot be reliably calculated because 
the physics processes responsible for bringing the neutralino relic density to its hnal value he 
outside the present framework. Put differently, we caution against unilaterally excluding model 
parameters (including the RNSb model) based on these considerations, because this frequently 
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Figure 6: Spin-in depen dent pZi scattering cross section vs. Mi (red dots) or M 2 (blue pluses) 
for the RNS benchmark point. 

requires other assumptions about the cosmological history of the Universe that have no impact 
upon collider physics]^ While WIMP discovery would be unambiguous, interpretation of the 
physics underlying any signal would require a careful specihcation of all underlying assumptions. 

The expected spin-dependent (SD) proton-neutralino direct detection cross section is plotted 
versus the gaugino mass parameter in Fig. In this case, the scattering occurs dominantly 
via Z-exchange. The ZZ^Zx coupling (Eq. 8.101 of Ref. |25]) is proportional to a difference 
in square of higgsino components of the neutralino. For Mi large and positive, both higgsino 
components are comparable and there is a large cancellation in the coupling. As Mi decreases, 
the higgsino components of Z\ decrease, but the up-type higgsino content more so than the 
down type. There is less cancellation and the coupling increases. As Mi decreases further, the 
bino component increases and the smallness of the higgsino components decreases the coupling. 
The negative Mi side shows more or less similar features until we reach Mi ~ —600 GeV where 
the flip in the identity of the neutralino mentioned in the previous hgure results in the jump. 

As far as WIMP detection goes, the SD cross section would influence IceCube ED detection 
rates the most since the WIMP abundance in the solar core is determined by equilibration 
between capture and annihilation of WIMPs in the sun, rather than on their local abundance. 
The scattering/ capture rate of the Sun depends mainly on the Hydrogen-WIMP scattering 
cross section which proceeds more through the SD interaction since there is no nuclear mass 

®What is clear from the data is that neutralinos with a large higgsino content (including the well-tempered 
neutralino) cannot be the bulk of the local dark matter. 
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Figure 7: Spin-dependent pZi scattering cross section vs. Mi (red circles) or M 2 (blue pluses) 
for the RNS benchmark point. 


enhancement. While some of the predicted values (red points) might well be marginally ex¬ 
cluded by the IceCube search, the take-away message is that for the most part the model with 
H = 200 GeV is on the edge of detectability, as long as neutralinos dominantly annihilate to W 
pairs. 

In Fig. we show the thermally-averaged neutralino annihilation cross section times rela¬ 
tive velocity evaluated as u —)■ 0. This quantity enters the halo WIMP annihilation rate, and 
detection rate for galactic positrons, anti-protons and gamma rays from WIMP halo annihila¬ 
tions are proportional to this factor. In the case of gaugino mass unihcation where we have a 
higgsino-like neutralino, then the local abundance is reduced and the expected detection rate 
is reduced by the square of the WIMP underabundance: where ^ /0.12. From the 

hgure, we see that while the local abundance increases as |Mi| is reduced (Fig. [^, the annihi¬ 
lation rate decreases because annihilation to WWs occurs mainly via the (reducing) higgsino 
component of the LSP. Once this channel is closed (around |Mi| ~ 200 GeV), annihilation to 
fermions takes over and the rate drops further. The FERMI-LAT collaboration has obtained 
upper limits located at about a few x 10“^®cm^/s (~ 2 x 10“^® cm?/s) for annihilation to bb 
{WW pairs) [5T]. Assuming a Navarro-Frenk-White prohle for dwarf galaxies in the analysis, 
models with a larger cross section would have led to a flux of gamma rays not detected by the 
experiment. Even without the scaling noted above, and certainly after the scaling, these 
bounds do not exclude any of the points in the hgure. For completeness we note that all the 
caveats that we discussed for the applicability of direct detection bounds are also applicable in 
this case, and we urge the reader to use caution in excluding ranges of parameters even if the 
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Figure 8 : Thermally-averaged neutralino annihilation cross section times velocity at n = 0 vs. 
Ml (red dots ) or M 2 (blue pluses) for the RNS benchmark point. 


Fermi Collaboration obtains tighter bounds in the future. 

3 Natural SUSY with a wino-like LSP 

In this section, we examine the phenomenological implications of altering the SU{2) gaugino 
mass parameter M 2 while keeping Mi = M 3 = 700 GeV. We begin by showing, as blue pluses, 
the variation of Aew with M 2 in Fig. Again, we see that Aew is relatively insensitive to 
M 2 except for the largest values of this parameter. This is due to the increasing contribution 
of winos to S“(hFi^ 2 )- Thus, models with M 2 Mi 3 lead to a wino-like LSP at little cost 
to naturalness. For M 2 < 150 GeV the chargino becomes lighter than 100 GeV (roughly the 
chargino mass bound from LEP2). Here, and in subsequent figures, we do not consider negative 
values of M 2 as these lead to a chargino LSP: 

In Fig. we show how the masses of charginos and neutralinos change as M 2 is reduced 
from its unified value. Starting with the RNSh spectra at M 2 = 700 GeV, where the W 2 and 
are essentially winos, and Zi, Z 2 and Wi are higgsinos, we see that as M 2 is lowered, the mass 
of the wino-like states reduces whereas the higgsino-like states remain with the mass hxed close 
to /r. The mass of the bino-like Z 3 also remains nearly constant. This behaviour persists until 
we reach the bino-wino level crossing near M 2 ~ 350 GeV where Z 3 and Z 4 switch identities. 
For still lower values of M 2 , we see another level crossing between the charged as well as neutral 
wino-like and higgsino-like states. For M 2 < 200 GeV, the lighter chargino as well as the LSP 
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Figure 9: Variation of chargino and neutralino masses vs. M 2 for the RNS SUSY benchmark 
model with variable M 2 but with Mi = M 3 


are wino-like, the heavier chargino and the neutralinos ^ 2,3 are higgsino-like, and Z 4 is mainly a 
bino. The mass gap ~'^Zi actually decreased with decreasing M 2 since these wino-like 
states have very tiny mass splittings. The mass gaps ^Z 2 ~^Zi g^^^citly increase 

with decreasing M 2 , reflecting the widening higgsino-wino mass difference. This should make 
their visible decay products harder so that these states are easier to detect at the LHC. 

We show the thermally-produced neutralino relic density versus M 2 in Fig. 
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Starting with M 2 = 700 GeV for which ~ 0.01, we see that steadily decreases 

with decreasing M 2 and reaches a value ~ 0.001 for very low values of M 2 where the 

Zi is nearly pure wino. This is because wino annihilation proceeds via the larger SU(2) triplet 
coupling to electroweak gauge bosons while annihilation of higgsinos proceeds via the smaller 
doublet coupling - the cross section for annihilation to W pairs, which is dominated by the 
f-channel chargino exchange, goes as the fourth power of this coupling. Thus, in the case of low 
M 2 with a wino-like neutralino, we might expect an even more reduced local abundance from 
thermally produced LSPs. The balance may be made up either by axions or other relics, or by 
LSPs produced by late decays of heavier particles. We cut the graph off when falls below 
its LEP2 bound. We do not see any dips corresponding to s-channel h or Z funnel annihilation 
as these fall in the LEP2 excluded region. 

An RNS benchmark point with a wino-like LSP is shown in Table [T] and is labelled as 
RNSw. All input parameters for RNSw are the same as for RNSh except now M 2 is chosen to 
be 175 GeV. The Wi — Zi mass gap has decreased to just 7.3 GeV while the Z 2 — Zi mass gap 
has increased beyond the RNSh value up to ~ 97 GeV, large enough so that both Z 2 ^ ZiZ 
and Z 2 —)■ decays are now allowed. In such a scenario, we would expect LHG SUSY 

cascade decay events to be rich in content of real Z bosons that could be searched for at the 
LHG. In fact, the GMS [52] and ATLAS [53] collaborations have already obtained bounds 
on chargino and neutralino masses from an analysis of about 20 fb“^ of LHG 8 data. These 
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Figure 10: Variation of vs. M 2 (blue curve) for the RNS SUSY benchmark model with 

variable M 2 but with Mi = M 3 . We cut the graph off at the low end because falls below 
its LEP2 bound. 

limits are obtained in simplihed models from an analysis of expectations from W 1 Z 2 and WiWi 
production at LHC 8 , assuming that and that the charginos (neutralinos) decay 

100% of the time to W bosons {Z bosons or Higgs bosons). While these limits are not directly 
applicable to pair produced W\ and Z 2 for the RNSw scenario in the table, the reader may be 
concerned that higgsino-pair production processes pp —)• W 2 Z 2 ,zX, W 2 W 2 X would lead to hnal 
states similar to what the LHC searches look for. We have checked that the RNSw scenario, 
with = 114 GeV is clearly allowed by current searches; data from the LHC13 run should, 
however, decisively probe this benchmark point. 


3.1 Implications for LHC13 

3.1.1 Gluino pair production: multijet plus events 


As discussed in Sec. 2 . 1 . 1 [ the discovery reach of LHC13 for gluino pairs mainly depends on the 
value of rUg which dictates the total gg production cross section in the case of heavy squarks. 
We would thus expect a similar LHC 13 reach for gluino pair production in the RNSw case 
as for RNSh and as for mSUGRA/CMSSM for comparable gluino masses and heavy squarks. 
Also, in the RNSw case, then charginos Wi will still be largely invisible due to their soft decay 
products. In some AMSB models with a wino-like LSP, then the mass gap lies at 

the 100 MeV level leading to long-lived winos whose tracks before decay may be visible 0 . 
In our case though , since fi is 100-200 GeV as required by naturalness, the Wi — Zi mass gap 
tends to lie in the 5-10 GeV range and so charged winos will be short-lived with no discernable 
tracks or kinks. However, in the RNSw case, then the Z 2 — Zx mass gap does become large and 
the well-known dilepton mass edge at should be observable for energetic enough 

Z 2 —>■ Zit^l~ decays if ~'^Zx ^ where the decay Z 2 —)• ZxZ opens up, then 
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the gluino cascade decay events (which, depending on the spectrum, should mostly proceed via 
real or virtual stop decays because stops are much lighter than hrst/second generation squarks) 
should be rich in OS/SF dileptons which reconstruct Mz- Note also that for modest values 
of M 2 , then Z 3 is also expected to be relatively light, and should also be accessible via gluino 
decays. For yet smaller values of M 2 , ^ 2,3 Zih may also be allowed and should occur with 
a comparable branching fraction to the decay to real Zs. 


3.1.2 Electroweak-inos at LHC 13 


In Fig. 11, we show NLO cross sections from Prospino [3l] for electroweak -ino pair production 
at LHC13 for the RNS benchmark but for variable M 2 . Chargino pair production - shown in 
the topmost frame - occurs via wino as well as via higgsino pair production. For large M 2 the 
latter dominates, but as M 2 is reduced, wino pair production increases in importance until it 
completely dominates for M 2 ~ 100 GeV. W 1 W 2 production, for the most part occurs via small 
gaugino/higgsino content, and so has a smaller cross section than the kinematically disfavoured 
W 2 W 2 production. The level crossing as the light chargino transitions from being higgsino-like 
to wino-like as M 2 reduces is also evident in the upper two curves. 

Chargino-neutralino production, shown in the middle frame, also occurs via wino as well as 
higgsino-pair production processes. For large values of M 2 , higgsino pair production dominates 
and WiZi ^2 production processes have the largest cross sections. For very small values of 
M 2 , pair production of winos is dynamically and kinematically favoured, and WiZi occurs at 
the highest rate. The higgsino-like states hF 2 ,^ 2,3 have masses /i and are also produced with 
substantial cross sections. Notice that W 1 Z 2 production remains signihcant even for small values 
of M 2 , presumably because it is favoured by kinematics (and increased parton luminosity). 

Neutralino pair production (shown in the bottom frame) can only occur via higgsino pair 
production since electroweak gauge invariance precludes a coupling of Z to neutral gauginos. 
As a result, Z 1 Z 2 production dominates for large M 2 . For small values of M 2 (where Zx 
becomes wino-like) Z 2 Z 3 production becomes important; however, Z 1 Z 2 production remains 
large because of large parton densities. 

We see that for M 2 ~ 300 GeV, the cross sections for W\Zx and W]W\ production processes 
increase rapidly with decreasing M 2 since W\ and Z\ become increasingly wino-like. However, 
since the W\ — Zx mass gap reduces even below the higgsino-LSP case, these states remain 
difficult - perhaps impossible - to detect. Possibly W^Wx production may be detectable via 
vector-boson fusion-like cuts in events where energetic jets with a large rapidity gap are required 
|55j . Although the cross section for wino-like WxZx production becomes very large at low M 2 , 
this process is difficult to detect. However, WxZ 2 production remains at viable rates even 
for low M 2 . In this case, one might look for relatively hard OS/SF dileptons from Z 2 decay 
recoiling against only soft tracks and E™®®. Other possibly more promising reactions at low 
M 2 include H 2 .^ 3 , W 2 Z 2 , ^’ 2^3 and maybe also Z 2 Z 4 production, since the decay products from 
both the chargino and neutralino should be relatively hard and can lead to events with 
three or more leptons, or real Z and Higgs bosons. As we mentioned, LHG collaborations are 
already searching for an excess of just such events [52l [53l |56]. Note also that WxZy, and ZxZ 2 
production each has a cross section in excess of 100 fb at low M 2 but would be considerably 
more difficult to detect. 
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Figure 11: Electroweak-ino pair production cross sections versus M 2 for the RNS SUSY bench¬ 
mark model with variable M 2 but with Mi = M 3 
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Figure 12 : Chargino and neutralino production cross sections with unpolarized electron and 
positron beams at a linear e’''e“ collider with ^/s = 500 GeV for the RNS SUSY benchmark 
model with variable M 2 but with Mi = M 3 


3.2 Implications for ILC 


At ILC, the natural SUSY scenario with low M 2 becomes both more challenging and richer. 


The cross sections for chargino and neutralino pair production at ILC500 are shown in Fig. 
for unpolarized beams. For M 2 = 700 GeV, we have the higgsino pair production reactions 
—>■ WiW{ and Z 1 Z 2 dominating. As M 2 is lowered, then the Wi becomes more wino-like 
and lighter leading to a larger cross section. However, the mass gap lUi — Zi drops below 
10 GeV making chargino pairs more difficult but likely still possible to detect with specially 
designed cuts. Beam polarization would serve to ascertain the higgsino/wino content of the 
chargino. Also, the Z 1 Z 2 reaction falls with decreasing M 2 as the Z — Z 1 — Z 2 coupling decreases 
iZ only couples to higgsino components). As M 2 falls below 300 GeV, the the Z 2 Z 2 , reaction 
turns on and grows in importance because the Z 3 becomes increasingly higgsino-like. Here, we 
expect Z 3 to decay via 2-body modes into Z-bosons or higgs bosons and Z 2 to decay either to 
2- or 3-body modes depending on the mass gap. This reaction should be distinctive and easily 
visible. 


3.3 Implications for WIMP detection 

WIMP detection for models with radiatively-driven naturalness and a wino-like WIMP may 
be either more and less difficult than the case with gaugino mass unihcation since, though 
the nucleon neutralino scattering cross section is larger, the local abundance for a thermally 
produced wino-like LSP is below the already low value typical of a higgsino-like LSP. Of course, 
the thermal wino abundance can be augmented by non-thermal processes involving moduli 
decay IS7I or axino/saxion decay m in the early universe. 

In Fig. we show the SI direct detection Zip scattering cross section versus M 2 as the 


20 







curve with blue pluses. Starting off at large M 2 , we see that as M 2 is decreased, the a^\Zip) 
cross section increases, and the increase is substantially larger than the case of a bino-like LSP. 
Recall this cross section proceeds mainly via light h exchange which depends on a product of 
gaugino and higgsino components of the neutralino LSP [25]. In this case, the wino-component, 
which involves the larger SU{2) gauge coupling g, becomes enhanced leading to the large cross 
section. For small enough M 2 < 250 GeV, the cross section turns around and decreases with 
decreasing M 2 since the Zi becomes more purely wino-like and the higgsino components are 
diminished. We note here that though the cross section in Fig. [^exceeds the stated bounds 
(1 — 2 X 10“® pb for = 100 — 200 GeV) in Ref. [29|, these bounds are not directly applicable 
because they were obtained assuming the neutralino constitutes the entire dark matter content 
of the Universe. For the natural SUSY scenario, the rates in direct detection experiments could 
be much smaller, as these scale by the neutralino fraction of the total local dark matter density. 
A wino-like neutralino that forms the bulk of the local dark matter would be excluded. 

In Fig. we show the spin-dependent direct detection cross section a^^{Zip) versus M 2 
as the blue curve. Here, the SD scattering cross section which proceeds mainly by Z exchange 
becomes large since there is less cancellation in the Z — higgsino — higgsino coupling. For 
small enough M 2 , then again the cross section turns over and decreases due to the diminishing 
higgsino components. We see that the cross section exceeds its 90% GL IceGube upper limit 
~ 1.5 X 10“^ pb [31] obtained assuming that LSPs in the sun annihilate dominantly to lU-pairs 
if M 2 < 700 GeV. Thus a truly wino-like neutralino scenario such as the RNSw benchmark 
is strongly disfavoured, as long as the assumption that the neutralino density in the sun is 
determined by equilibrium between capture and annihilation rates is valid. 

In Fig. I^we show {av) versus M 2 as the blue shaded curve. In this case, as M 2 falls, then 
Z\Z\ —)■ WW becomes large and the annihilation rate increases. One might expect increased 
liklihood for indirect WIMP detection via gamma rays and antimatter detection. However, the 
increased annihilation rate is counter-balanced by a likely decreasing local WIMP abundance 
where the detection rate is proportional to the square of the reduced local abundance. We see 
that although the predicted rate naively exceeds the upper limit from Fermi-LAT in Ref. EH, 
after the scaling discussed above exclusion is not possible. 

4 General results in Mi vs. M 2 plane 

While it is instructive to examine natural SUSY models with reduced GUT scale bino- or wino- 
mass parameters, there is no compelling reason to believe that one parameter is unified with 
M 3 (GUT) while the other is quite different. In general one may have arbitrary gaugino masses 
and in fact both may be reduced leading to a mixed bino-wino-higgsino LSP. Here, we present 
some illustrative studies of this more general situation. We can choose Mi > 0 by convention. 
The signs of M 2 and M 3 as well as p are then physically relevant. Since our purpose is to give a 
broad brush idea of how RNS phenomenology of electroweak-inos may be altered, we will take 
M 3 and p to be hxed at their values for the RNSh benchmark point and display results in the 
Ml — M 2 planep^ 

In Fig. we show the Mi vs. M 2 plane for the RNS benchmark model but with Mi and 
^°The electroweak sector should be almost insensitive to M 3 , but will, of course, be sensitive to the sign of g. 
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Figure 13: Dominant component of the neutralino LSP in the Mi vs. M 2 plane for the RNS 
SUSY benchmark model. The LSP is dominantly a bino, wino or higgsino in the region denoted 
by blue dots, red pluses and green crosses, respectively. Other parameters are hxed at their 
values for the RNSh model point in Table [T| In the region marked LEP2 excluded, < 
100 GeV, whereas in the remaining unshadea region of the lower half plane, 


M 2 as free parameters. The black dot in the upper-right corner denotes the location for unihed 
gaugino masses. The regions of the plot are coded according to the dominant content of the 
Zi. bino (blue dots), wino (green triangles) and higgsino (red pluses). The special cases of 
the previous sections correspond to moving horizontally to the left or vertically down from the 
unihed gaugino mass point. We start the scans at Mi = 50 GeV and scan both signs of M 2 . 
Here, and in subsequent hgures, the band with IM 2 I ~ 150 GeV is excluded by the LEP2 bound 
on the chargino. In the half plane with M 2 < 0, the additional region without any shading 
corresponds to a charged LSP (m^^ < and so is excluded by cosmological considerations. 
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we show the Wi — Zi mass gap in the Mi vs. M 2 plane for the RNS benchmark 


In Fig. 

model. The purple shaded region has mass gaps between 10 and 20 GeV and corresponds to 
the bulk of the higgsino-like LSP region along with the wino-like LSP region. As expected, the 
mass gap becomes small when lUi and Zi are both higgsino-like (|p| ^ |Mi^ 2 |) or when these 
are both very wino-like (M 2 <C |p|). It also becomes small along the boundary of the region 
in the lower half plane where the chargino becomes the LSP. It is mainly when one moves to 
small Ml, or large IM 2 I and moderate Mi, that this mass gap exceeds 40-50 GeV, so that the 
daughter leptons from chargino decays are expected to be relatively hard. 

The Z 2 — Zi mass gap is shown in the Mi vs. M 2 plane in Fig. 15 Typically the smallest 


mass gap occurs when we have a higgsino-like LSP as in the case of gaugino mass unihcation 
in the upper right part of the plane, or in the region where Mi and IM 2 I are both much larger 
than /i (M 2 < 0). A small (purple) mass region also occurs when Mi ~ I 2 M 2 I (lower half plane) 
so that the weak scale values of Mi and IM 2 I become comparable upon renormalization group 
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Figure 14: The mass gap in the Mi vs. M 2 plane for the RNS SUSY benchmark 

model. 


^1 


evolution: i.e. the bino and wino become nearly degenerate, but the states remain nearly pure 
winos and binos because of opposite signs of their mass terms. In this very low Z 2 — Zi mass gap 
region one might expect enhanced bino-wino co-annihilation(BWCA) in the early universe [58] . 
Note that the Z 2 — Zi mass gap in especially the upper half plane, exceeds 50 GeV for a large 
swath of the plane, and is larger than Mz and even nih over a substantial part. This should 
make for interesting signals at LHC13 via the multilepton, WZ and Wh plus channels at 
LHC 13. It is also noteworthy that a region exists where both the lUi — Zi and Z 2 — Zi mass 
gaps fall below 10 GeV. This occurs in the narrow crescent at large Mi in the lower half plane. 
This region might be challenging even at the ILC if the heavier charginos and neutralinos are 
kinematically inaccessible. In this case, techniques using initial state photon radiation might 
be required [59] . 

We note that while we have focussed on the mass gap between the lighter charginos and 
neutralinos and the Zi, there is a substantial region of the parameter space of natural SUSY 
models where signals from the heavier charginos and neutralinos should be accessible at LHC13. 
Because CMS and ATLAS LHC searches [52l [53] tend to employ hard cuts, it is entirely 
possible that signals from the heavy states (assuming these are within the LHC13 reach) reveal 
themselves more easily than signals for the lighter states. 

The thermally-produced neutralino relic density is shown in Fig. 


16 The regions with very 


low relic density ~ 0.01 are 1. the wino-like LSP region along with 2. the BWCA strip 

in the lower half plane where — ^z^ resonance annihilation regions where 

~ Mz or rrih (the vertical strips at low Mi). The thermal relic density is also below 
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Figure 15: The mass gap in the Mi vs. M 2 plane for the RNS SUSY benchmark 

model. 

its observed value in the higgsino region or in parts of the mixed bino-higgino LSP region. In 
these regions, we will need either additional dark matter particles or non-thermal production of 
neutralinos to match the measured value of cold dark matter relic density. The boundary of the 
light- and dark-blue shaded region is where we have a well-tempered neutralino whose thermal 
neutralino relic density can saturate the cold dark matter. In the light-blue and green-shaded 
parts of the plane (deep in the bino LSP and away from the Z and h resonances) the relic 
density of neutralinos must be diluted by entropy production late in the history of the Universe 
or else the Zi must be made to decay either via i?-parity violating interactions or decay to an 
alternative LSP {e.g. an axino). 

We do not show the dark matter detection cross sections in this plane, partly because for the 
most part we do not expect that these will unambiguously constrain the parameter regions for 
reasons that we discussed earlier regarding the assumed local density of WIMPs. A noteable 
exception may be any constraints on the neutralino from the spin-dependent cross-section 
obtained by IceCube as these are valid as long as there is equilibration between neutralino 
annihilation rates and gravitational capture of WIMPs by the sun. 

5 Conclusions: 

Supersymmetric models with radiatively driven naturalness are especially interesting since they 
allow for Mzi rrih ~ 100 GeV whilst sparticles other than higgsinos can naturally be at the 
multi-TeV scale. Such spectra seem to be required by reconciling naturalness with LHC8 
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Figure 16: Thermally-produced neutralino relic abundance in the Mi vs. M 2 plane for the RNS 
SUSY benchmark model. 

sparticle search constraints and with the measured value of the Higgs boson mass [001 El] • Most 
previous analyses have examined RNS models in the context of gaugino mass unification. In that 
case, the LSP is expected to be higgsino-like and constitute only a portion of the dark matter 
while axions could make up the remainder. The light higgsinos required by naturalness can 
evade LHC searches because of their compressed spectrum: higgsino decays release only small 
visible energy, so that their production remains hidden under Standard Model backgrounds. 

These results follow from requiring both naturalness and gaugino mass unihcation. We regard 
naturalness to be one of the main motivations for supersymmetry. In contrast, while gaugino 
mass unification is highly motivated by the simplest GUT models, it is easy to construct GUTs 
with non-universal gaugino masses at no cost to naturalness. Gaugino mass non-universality 
results if vacuum expectation values of the auxiliary fields that spontaneously break super- 
symmetry also break the GUT symmetry. The main requirement from LHG searches is that 
M 3 ~ rrig ~ 1.3 TeV. The values of bino and wino mass parameters are relatively unconstrained. 
If their weak scale values are similar to or less than |/r|, then the LSP can be either bino-like or 
wino-like (or a mixture) instead of just higgsino-like at no cost to naturalness. In such a case, 
both the collider expectations and dark matter/WIMP search expectations change in important 
ways. 

We have shown that in the case of natural SUSY models with enhanced bino LSP content, 
increased mass gaps Wi — Zi and Z 2 — Zi are expected on account of bino-higgsino mixing. 
The harder decay products of Wi and Z 2 lead to discernable effects such as the presence of 
dilepton mass edges in LHG events, and perhaps additional light electroweak -ino pair produc- 
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tion processes at the ILC involving also Z 3 production. In the wino-like LSP case, then only 
the Z 2 — Z\ mass gap opens up, while Wi — Z\ gap becomes tighter. This situation should 
be readily discernable at ILC, especially with the availability of polarized beams. Of course, 
if M 2 and |/i| both assume modest values, the heavier states IV 2 , ^ 3,4 will also be accessible 
at the LHC, and electroweak chargino and neutralino production will lead to a rich variety of 
multilepton, WZ and Wh plus events that are already being searched for [521 ESI ES], 

and possibly also spectacular events without additional jet activity. In such a 

scenario, ILC would become both a higgsino and a wino/bino factory and it should be possible 
to perform a detailed bottom-up study of the electroweak-ino sector, assuming that all states 
are kinematically accessible [62] . 

Expectations for WIMP searches also change. In the case of a bino-like LSP, we generally 
expect a larger thermal abundance of neutralino dark matter. While it is possible to obtain a 
well-tempered neutralino that saturates the observed cold dark matter relic density, the thermal 
neutralino density is often too large in which case it needs to be diluted by late-time entropy 
production or else allowed to decay. As a result, the neutralino contribution to the relic density 
today depends on the (unknown) physics, leading to signihcant uncertainties in prediction of 
rates for direct detection searches. While this makes it difficult to use experimental bounds from 
LUX/XENON100[29| and other experiments to unambiguously exclude portions of parameter 
space without a complete model of particle physics and cosmology, these searches could lead 
to a discovery! 

In the case of natural SUSY with a wino-like WIMP, then one expects an even lower local 
abundance from thermally-produced neutralinos as compared to the value for higgsino-like 
LSPs. The measured relic density must then be made up by other (non-WIMP) relics of 
which axions may be the most promising, or via WIMP production from late decays of heavy 
particles. In view of the resulting uncertainty in the expectation for local density of neutralino 
dark matter, we once again advocate using caution when interpreting the absence of events in 
direct and indirect dark matter searches to exclude ranges of model parameters. The IceCube 
upper limits inn on the spin-dependent neutralino-nucleon scattering cross section may be 
an exception in this regard, because the validity of these bounds depend on the assumed 
equilibration between the neutralino capture and annihilation rates in the sun, and not on the 
local density of dark matter. Indeed portions of parameters with low values of M 2 may already 
be excluded since we expect that equilibration has already occurred. 

To sum up, in our view, supersymmetric GUTs remain the most attractive solution to 
the naturalness problem plaguing the Standard Model and light higgsinos are the most robust 
consequence of naturalness considerations. If electroweak gaugino mass parameters happen to 
assume modest values - this is not required by naturalness but is completely compatible with 
it - there could be spectacular signals from electroweak gaugino production at the LHC in 
multi-lepton-|-E™^®®, lUZ-l-E™®®, IUh-|-E™®® and -t-E™'®® channels. Direct and indirect 

searches for WIMPs could also reveal a signal even in the case of a depleted local abundance 
of WIMPs. If natural supersymmetry is realized with fortuituously low gaugino masses, then 
prospects for SUSY discovery at LHC13 will be vastly improved since signals from several 
chargino and/or neutralino states might also be observable. Production of light electoweak -ino 
states at ILC - as required by naturalness [H] - remains true but with even richer prospects 
since both gauginos and higgsinos could be kinematically accessible. 
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